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Abstract 

We present a lattice investigation of heavy-light mesons in the quenched 
approximation, using non-relativistic QCD for the heavy quark and a clover 
improved Wilson formulation for the light quark. A comprehensive calcu- 
lation of the heavy-light spectrum has been performed for various heavy 
quark masses around the b. Our results for the Bg — splitting agree well 
with the experimental value. We find the Aj, — B splitting to be compat- 
ible with experiment, albeit with large error bars. Our B* — B splitting 
is slightly low, which could be explained as an effect of quenching. For 
the first time, we are able to estimate the mass of P states at the B and 
compare them with experiment. 
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1 Introduction 



Heavy-light hadrons containing a b quark are of topical interest in particle theory. 
Theoretical predictions for heavy-light meson decay constants are needed for an 
experimental determination of elements of the CKM matrix A theoretically 
interesting feature of heavy-light hadrons is that, in the limit of a very heavy 
quark mass mg, the dynamics simplify and can be described by an effective field 
theory (Heavy Quark Effective Theory) P|. In the heavy quark limit mq — > oo 
additional symmetries between states of different spin and flavour appear which 
make it possible to relate different decay processes to each other and to make 
qualitative predictions about the spectrum. Away from the heavy quark limit a 
systematic expansion in the heavy quark mass can be performed whose coefficients 
have to be determined non-perturbatively. On the lattice this can be done from 
first principles. 

Since, for the values of the lattice spacings which can be achieved in present 
computer simulations, the b quark mass is > 1 in lattice units, b quarks cannot be 
simulated directly using a naive relativistic quark action. It has been suggested 
that heavy quarks around and above the charm can be simulated using a re- 
formulation of a relativistic quark action Q], but there are indications that there 
are problems with this method in the b quark region at least in the way 
the method has been implemented to date. So simulations are often done at 
smaller quark masses, up to the charm, and results extrapolated to the b. With 
this method it is very difficult to get certain spectral quantities, such as spin 
splittings, to agree with experiment. Simulations of heavy-light systems in the 
static (i.e. infinite mass) limit of the heavy quark give sensible results only 
for spectral quantities that are fairly independent of mq. For other quantities, 
contributions due to the finiteness of the heavy quark mass have to be calculated 
directly as perturbations in l/mq (see e.g. |0]). In addition the signal-to-noise 
ratio is much worse than for propagating quarks. 

The non-relativistic formulation of QCD (NRQCD) p however allows us, in 
principle, to simulate at the b quark mass. The largest momentum scale which 
governs discretization errors is not the heavy quark mass but the non-relativistic 
momentum of the heavy quark, which is much smaller. This method has been 
successful in precision spectroscopy of heavy- heavy systems such as T and J/ "if 
mesons §]. In these systems the energy scale inside the meson is set by the 
kinetic energy of the heavy quarks, which is of the order of Aqcd- Relativistic 
corrections are included in a systematic expansion in Vq, where vq is the velocity 
of the heavy quark. For mesons with one heavy and one light quark the power 
counting rules for the non-relativistic expansion are determined by the fact that 
here the energy scale for the heavy quark inside the meson is the binding energy. 
For rriQ oo only the time derivative appears in the tree level Langrangian. 
The covariant time derivative acting on the heavy quark spinor gives the quark 
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energy: 

DtQ ~ ^bindQ, (1) 

Due to momentum conservation the heavy and hght quark momenta are equal in 
the rest frame of the meson: 

mgv = pQ=pgr^ 0{Aqcd), (2) 

where Aqcd takes a value around 300 — 500 MeV. So the heavy quark velocity 
is small: 

v,-0 (^] , (3) 

which is of the order of 10% for the B, and we can choose vq, or alternatively 
l/mq, as the expansion parameter for the NRQCD interactions. Now we have 
to determine the order in v of various correction terms. We can immediately 
write down the operators which appear with a tree level I/ttlq coefficient: the 
non-relativistic kinetic energy operator 

Okin = (4) 

2mQ 

and the magnetic interaction operator 

To estimate the size of the heavy quark kinetic energy, we write 

0^ ^ {mgyf ^ Aqcd .gx 
2mQ 2mQ 2mQ ' 

thus we expect Okin to be suppressed with respect to Dt by a factor of Aqcd/'^i^q- 
For the magnetic operator we have 

5 = V X A (7) 

Since 

A~p~0(Aqc7d) and V ~ p ~ O(Aqcd), (8) 
we estimate the magnetic operator to also be of order 

~ -2^- 

The heavy quark formulation used in this calculation is described in more 
detail in section where we also talk briefly about the implementation we use for 
the light quark and discuss the choice of our simulation parameters. In section ^ 
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we show the operators from which our meson and baryon states are constructed. 
Details of our fitting procedure and the results are discussed in section |. We 
also describe there how meson masses are extracted and give results for the Bd 
and the Bs — Bd splitting. After that we turn to the B* — B (hyperfine) and 
Ab — B splittings and compare simulation results with experimental values and 
expectations from HQET. Finally, we present preliminary results for P wave 
states. The conclusions are given in section |^. 



2 Simulation details 

The gauge fields used here are an ensemble of 35 quenched configurations at 
P = 6.0, generated by the UKQCD collaboration. The configurations are sep- 
arated by 200 sweeps, each including 5 overrelaxation and 1 Cabibbo-Marinari 
heatbath sweep. By time reversing them, we obtain an additional set of 35 con- 
figurations which we treat in our analysis as statistically independent. We fixed 
the configurations to Coulomb gauge. 

In this simulation the non-relativistic Lagrangian, which describes the b quark, 
is expanded through order 1 / mg at tree level: 

C = Q^iDt + Ho + 5H)Q, (10) 

where 

H,^-^,5H = -S^ (11) 

To get rid of large renormalizations of the link operator, we use tadpole improve- 
ment |T^ : 

^ U^/uo, 4 = {^.TrUpi,^.) , uo - 0.878, (12) 

so we can use the tree level coefficient cb = 1- The heavy quark propagator is 
computed using the following evolution equation 

and on the following timeslices 

n / „ TJ \ n 



= [l- ^yul (l - (1 - a6H)G„ t > 1, (14) 

for the heavy quark masses amq = 1.71, 2.0, 2.5 and 4.0. The bare heavy quark 
mass corresponding to the b quark should be the same as the one tuned in T 
spectroscopy calculations with NRQCD. If the b quark mass is set to am° = 1.71, 
the simulated T mass agrees with the experimental value [0], using a lattice 
spacing of 2.4(1) GeV, obtained from the IP — IS and 2S — IS" splittings of the 
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T. However, due to the effect of quenching one obtains different lattice spacings 
from different physical quantities, or in other words, for different physical systems 
different scales are appropriate. From light hadron spectroscopy at the same value 
of P one gets lattice spacings closer to ~ 2 GeV. Probably, scales in heavy- 
light systems are closer to the ones in light hadrons. In this paper we use 2.05(10) 
GeV. The central value is taken from the APE collaboration who performed a 
spectrum calculation at the same value of jS and the same type of light fermions 
(clover fermions with unity clover coefficient) and quote a value of = 2.05(6) 
GeV from rrip P]. We use a larger error for a~^, to encompass other results from 



^ 2 from nip 



light spectroscopy (e.g. a = 2.11(11) from /^r in 0, a ^ = 2.0 
in flOl-)) well as from the string tension |jll| of = 1.94(8) GeV at /? = 6.0. 
Systematic errors in light hadron spectroscopy are still rather large and not well 
determined, so the uncertainty in may be underestimated. For a^^ = 2.05 
GeV, arriQ = 2.0 corresponds to the bare b quark mass used in the T simulations. 
The stability coefficient n has been set to 2 for all masses. 

We also compute the heavy quark propagator in the zeroth order (static) 
approximation of the l/mq expansion, where the Lagrangian is just given by the 
covariant time derivative: 

C = Q^DtQ. (15) 
and the heavy quark propagator follows the evolution equation: 



G 



UlG, 



5, 



{x,t)fl- 



(16) 



For the light quarks we use propagators generated by the UKQCD collabora- 



tion with a clover action 12 



E 



+ + /i)(l + -fp)Ul{x)q{x] 



H — q{x)q{x) 



lean 



/2j2<l{x)(T^uFf,^q{x) 



(17) 



with clover coefficient c = 1, i.e. without tadpole improvement. The quark fields 

are rotated at the source: 



q{x) 



1 — ■ D ] q{x 



It has been shown that matrix elements calculated with these quark fields 
are free of 0(a) errors at tree level. 

We use light fermions at values of k, of 0.1440 and 0.1432, which bracket the 
strange quark mass. On the same set of configurations and light propagators the 
critical k values has been determined to be 0.14556(6). The strange quark, whose 



mass is determined from mj^/rn^, corresponds to Kg = 0.1437 1 5 [TDI. 
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3 Hadron operators 



Our meson states are built up from various interpolating operators: 

Y.Q\xi)T\xi-X2)q{x2). (19) 



Xx,X2 



r(xi — X2) factorizes into a smearing function (/)(r = \xi — X2\) and an operator 
VL which consists of a 4 x 2 matrix in spinor space and, for P states, a derivative 
acting on (p. We choose to be either a delta function or, to project onto 
the ground or an excited state, a hydrogen-like wave function. The smearing is 
applied on the heavy quark. In general the smearing function is different at 
the source (sc) and at the sink (s/c), and we obtain meson propagators of the 
following form: 

C{p = 0,t)=Y: Tr [75(M-i)t(y2)75r(^^')t(yi -y2)G,(yi)] • (20) 
yi,y2 

Here, is the light quark propagator and 

Gi(f)=5:G,(y-f)r(-)(f) (21) 

X 

is the heavy quark propagator smeared at the source. 

The continuum quantum numbers of the S and P states we implemented and 
the corresponding lattice operators Q can be found in table |I|. For the 7 matrices 
we chose the Dirac basis using the representation 





f 1 






f 


—ia \ 


7o = 


K 


-J 




I ia 


j 



(22) 

Note that in heavy-light systems C is not a good quantum number, and we expect 
mixing between the ^Pi and the ^Pi states. 

For the baryon we used the following interpolating operator: 

Oa = eabcT.Qai^i)T.iiyi^2)Cj5loq'{S2mx, " ^2!), (23) 

XI X2 

where C = 7072 is the charge conjugation matrix. The baryon correlators are 
smeared at the source and local at the sink. 



4 Results 

In the following we will use the notation Crs for the correlation functions, where 
the index r denotes a delta function (L), ground state smearing function (1) 
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Meson Lattice 
(JP) Rep. 



n 



'So (0-) Ai 



'S. (1-) 
ip, (1+) 
'Po (0+) 
'Pi (1+) 
=^^2 (2+) 



1» 



^i(i) 



C) 

(i) 







(i ^ 3) 



Table 1: Meson Operators. A j denotes the symmetric lattice derivative. 1 stands 
for the 2x2 unit matrix. 
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or excited state smearing function (2) at the source. The index s denotes the 
smearing function at the sink. 

From the original data we generated 100 bootstrap ensembles, each containing 
the original number of samples, and fitted to each of them separately. The fit 
results and errors are obtained from averaging over the number of bootstrap 
samples. This procedure also enables us to take correlations between data with 
different k and rriQ values from the same configuration into account. Given the 
relatively low statistics, we are sometimes forced to discard eigenvalues in our 
SVD inversion of the covariance matrix. This is done using a cutoff A on the 
ratio between the smallest and the largest eigenvalue. 

4.1 Ground state energy 

To extract the bare ground state energy, -Esim, and amplitudes we fitted the 
correlators Cn and Cn simultaneously to a single exponential. The signal is good 
up to large times and we fit out until t^ax = 25. The plateau of the correlation 
functions at k = 0.1432 is reached at t^m = 5, but at k = 0.1440 there is a 
slight decrease in the ground state energy and especially the amplitudes until 
train is movcd out to 9. Moving tmin further out does not change the results 
any more. An example for the dependence of the ground state energy on the fit 
interval at the different n values is shown in table 0, together with the quality of 
fit parameter Q. 

With our limited statistics, our bootstrap procedure generates certain en- 
sembles on which multi-exponential fits fail, but with the original ensemble of 
correlation functions we can do a simultaneous fit of Cil and C2L to 2 exponen- 
tials. This gives a value for the ground state energy which tends to be a little 
higher, but in general still compatible within one standard deviation with the 
results fitted with just one exponential. Results of double exponential fits for 
arriQ = 4.0 are shown in table |^. The single exponential fits give higher Q values 
and better looking effective mass plots, so we use them to extract the meson 
ground state energies. We can also make a rough estimate of the energy of the 
first excited state. Its error bars are too large to extract a dependence on the 
heavy or the light quark masses. From the results in table ^ and fit results at 
other heavy quark masses we take the value of £'sim(2S') = 0.8(1) as a reasonable 
estimate. For a more reliable calculation of the excited state energy one would 
have to do three exponential fits to Cn and C2L, but, given our low statistics, 
errors on these fits are very large. 

We extract meson energies from fits in the range tmm/tmax = 9/25 for both 
K, values. The systematic error due to the variation of the fit interval and the 
difference between the single and double exponential fits is smaller than the fit 
error at this fit range. Effective mass plots for Cn and Cu are shown in figure |T[ 
Results for the pseudoscalar, vector and spin-averaged ground state energy E = 
(3-E'sim(^5'i) + i?sim(^'S'o))/4 for all heavy and light quark mass values are shown 
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K = 0.1432 


K = 0.1440 


tmin/tmax 


aEsim(15') 


Q 


a^sim(15') 


Q 


5/25 


0.516(6) 


0.3 


0.497(6) 


0.2 


6/25 


0.515(7) 


0.4 


0.495(7) 


0.2 


7/25 


0.514(6) 


0.4 


0.493(8) 


0.2 


8/25 


0.514(7) 


0.4 


0.491(7) 


0.3 


9/25 


0.513(7) 


0.4 


0.491(8) 


0.3 


10/25 


0.511(5) 


0.4 


0.489(7) 


0.3 



Table 2: Results of single exponential fits to the ground state correlation functions 
Cii and CiL at am% = 4.0. For k = 0.1432, A = 0.03, for k = 0.1440, A = 0.02. 







= 0.1432 






= 0.1440 






aE^i^{lS) 


a^sim(25) 


Q 


aEsim(l'S') 


a^sim(25) 


Q 


4/24 


0.513(4) 


0.81(4) 


0.07 


0.495(6) 


0.79(5) 


0.09 


5/25 


0.518(5) 


0.77(6) 


0.17 


0.497(7) 


0.71(6) 


0.19 


6/25 


0.516(5) 


0.79(8) 


0.15 


0.495(9) 


0.71(9) 


0.15 


7/25 


0.519(3) 


0.94(34) 


0.09 


0.502(5) 


0.92(25) 


0.08 


8/25 


0.520(7) 


0.91(18) 


0.07 


0.498(3) 


0.92(43) 


0.11 


9/25 


0.517(7) 


0.90(18) 


0.06 


0.501(5) 


0.88(25) 


0.09 


10/25 


0.517(6) 


0.88(16) 


0.03 


0.501(2) 


0.86(8) 


0.07 



Table 3: Results of double exponential fits to the correlation functions Cn and 
C2L for am° = 4.0, A = 10"^ 



in table ^, the chiral extrapolation and interpolation to the strange light quark 
mass in table H. 

The dependence of the bare spin-averaged ground state energy on the bare 
heavy quark mass is illustrated in figure 0. This quantity is spin-independent and 
thus useful to estimate the correction due to the kinetic operator Okin- Correlated 
fits of the simulation results to linear functions give a slope of —0.02(2) in lattice 
units for k = 0.1440 and for k = 0.1432, —0.02(1). In both cases, the extrap- 
olation to the static limit is compatible with the static simulation results. We 
note that in the tadpole improved case, Okm is not a positive definite operator. A 
more physical quantity than Egim is the binding energy, -Ebind = -^sim — Eq, where 
Eq is the shift in the zero of energy for NRQCD or the static theory, respectively 
(perturbative values from []T3[ listed in table This quantity rises slightly as 
a function of I/ttiq, but there are significant perturbative uncertainties. For a 
more detailed discussion see 1116 . 
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«; = 0.1432 


K = 0.1440 




w 


n E ■ 


Q 


n E ■ 


Q 


'So 


1.71 


0.501(9) 


0.3 


0.480(7) 


0.4 




2.0 


0.506(6) 


0.4 


0.483(7) 


0.4 




2.5 


0.511(5) 


0.4 


0.487(7) 


0.4 




4.0 


0.513(7) 


0.3 


0.491(8) 


0.3 




oo 


0.524(6) 


0.2 


0.508(7) 


0.2 




1.71 


0.521(6) 


0.3 


0.502(7) 


0.4 




2.0 


0.523(5) 


0.3 


0.503(6) 


0.3 




2.5 


0.524(6) 


0.3 


0.504(6) 


0.3 




4.0 


0.527(6) 


0.3 


0.507(8) 


0.2 


spin avg. 


1.71 


0.516(6) 


0.3 


0.497(7) 






2.0 


0.518(5) 


0.3 


0.498(7) 






2.5 


0.521(6) 


0.3 


0.499(7) 






4.0 


0.524(6) 


0.3 


0.503(7) 





Table 4: Fit results for the bare ground state energy for the pseudoscalar and the 
vector meson and the spin-averaged energy for both k. values. 





ruq = 


TTLq ''^strange 




aniQ 




n E ■ 


'So 


1.71 


0.437(19) 


0.489(6) 




2.0 


0.437(11) 


0.491(6) 




2.5 


0.439(12) 


0.496(7) 




4.0 


0.446(15) 


0.499(6) 




oo 


0.461(15) 


0.511(6) 


'S^ 


1.71 


0.464(15) 


0.509(6) 




2.0 


0.463(12) 


0.510(6) 




2.5 


0.464(9) 


0.512(6) 




4.0 


0.478(12) 


0.515(6) 


spin avg. 


1.71 


0.457(14) 


0.504(6) 




2.0 


0.456(12) 


0.505(6) 




2.5 


0.457(9) 


0.508(6) 




4.0 


0.462(13) 


0.511(6) 



Table 5: Bare ground state energy, extrapolated to a vanishing (left) and the 
strange (right) light quark mass. 
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Figure 1: Effective masses of Cil (left) and Cn (right) correlation functions at 
am% = 2.5 and k = 0.1432. 
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Figure 2: Bare spin-averaged ground state energy as a function of the inverse 
bare heavy quark mass amg. Squares denote k. = 0.1432 and circles k, = 0.1440. 
The hnes represent correlated fits of the NRQCD results to a linear function. 
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n 


aEo 


aA 


aMp 


MpJGeV] 


K = 0.1432 


K = 0.1440 


TTlq = 


1.71 


2 


0.32(6) 


1.73(13) 


2.23 


2.21 


2.17 


4.5(5) 


2.0 


2 


0.30(6) 


2.02(12) 


2.53 


2.50 


2.46 


5.1(5) 


2.5 


2 


0.29(6) 


2.48(9) 


2.99 


2.97 


2.92 


6.0(5) 


4.0 


2 




4.07(5) 


4.58 


4.56 


4.52 


9.3(6) 



Table 6: Energy shifts used in our calculation and pseudoscalar meson masses 
calculated with these shifts. The last column contains meson masses for rUg — > 
in physical units, using a lattice spacing of = 2.05(10) GeV. Errors on aMp 
are dominated by the perturbative error in A, errors on the physical masses by 
perturbation theory and the uncertainty in the lattice spacing. 



4.2 Meson masses 

We calculate the meson mass using the relation 

Mb = A + (24) 

The energy shift A can be determined in perturbation theory. It consists of the 
renormalized heavy quark mass and a shift in the zero of energy: 

A = Z„m° - Eo. (25) 

The quantities Zm and Eq are calculated at 0{a): 

aEo = Aav{q*A), Z^ = I + Bav{q*B), (26) 

where we take for ay the two loop expression for the running coupling constant 
in the Lepage-Mackenzie scheme. The one loop terms A and B and the appro- 



priate momentum scales q\ and were calculated by C. Morningstar The 
shifts used in our calculation are shown in table together with results for the 
meson masses in lattice and physical units. Eq values are perturbative, shifts are 
perturbative except in the last row. At = 4.0, a non-perturbative energy 
shift calculated from the T has been used, to avoid defects in the scale q^ in 
the heavy quark mass renormalization, occuring at heavy quark masses around 
arriQ = 5. It is possible to calculate the meson mass non-perturbatively from 
the meson dispersion relation [|l7l. The non-relativistic dispersion relation reads 
in lowest order: 

Eip) = Eip = 6) + ^. (27) 

Note that E{p = 0) is, by definition, -Egim- We calculate the 'kinetic' meson mass 
Mb from the energy splitting between mesons of momentum \p\ = 2tt/Ls and 
IpI = 0, where Ls is the spatial lattice extent. 
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K = 0.1432 


K = 0.1440 




aAE Q aMp 


aAE Q aMp 


1.71 
2.0 
2.5 
4.0 


0.030(7) 0.5 2.6(6) 
0.030(12) 0.5 2.6(1.0) 
0.024(8) 0.5 3.2(1.1) 
0.018(8) 0.4 4.3(1.9) 


0.034(5) 0.5 2.3(3) 
0.030(5) 0.5 2.6(4) 
0.026(5) 0.5 3.0(6) 
0.018(5) 0.5 4.3(1.2) 



Table 7: Splitting between correlation functions with \p\ = 2ti/Ls and |p| = 
and meson masses derived from it. 

In our analysis of the energy splittings between the finite and zero momentum 
correlation functions we make use of the fact that correlation functions of different 
states on the same configuration are highly correlated. We generate 100 bootstrap 
measurements of the jackknifed ratio R of the correlation functions of the two 
states and fit these to a single exponential: 

R ~ Ae"^^*, (28) 

where AE is the energy splitting. 

As shown in table the masses obtained with this method have much larger 
errors, but are compatible with the masses calculated using perturbation theory. 

4.3 Mass splittings 
4.3.1 B, - Bd splitting 





aAE 


1.71 


0.047(11) 


2.0 


0.049(8) 


2.5 


0.050(5) 


4.0 


0.049(7) 


oo 


0.051(13) 



Table 8: Spin-averaged Bg — Bd splitting. 

Results for the spin-averaged B^ — B^ splitting are listed in table ||. As 
expected from HQET, this splitting is fairly independent of the heavy quark 
mass. Figure ^ shows the splitting as a function of the inverse spin-averaged 
meson mass with a correlated fit to a constant. The fitted value for this constant 
is 0.049(5), with a P^i' degree of freedom of 0.03. Converted into physical 
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0.07 



0.06 

m 

m" 0.05 



0.04 





1 


[ 


] 








: 


] 


c 


] 

[ 


] 



0.03 ^ — ' — ' — ' — ' — ^ — ' — ' — ' — ' — ^ — ' — ^ 
0.0 0.2 0.4 

l/(aM) 

Figure 3: Spin- averaged energy splitting between Bg and Bd in lattice units 
plotted against the inverse spin-averaged meson mass. The dotted line denotes a 
correlated fit of the simulation results to a constant. 





K = 0.1432 


K = 0.1440 




aAE 


Q 


aAE 


Q 


1.71 


0.0201(17) 


0.43 


0.0200(24) 


0.42 


2.0 


0.0180(14) 


0.51 


0.0180(20) 


0.47 


2.5 


0.0151(12) 


0.43 


0.0150(17) 


0.37 


4.0 


0.0103(10) 


0.28 


0.0102(14) 


0.22 



Table 9: Fit results for the B* — B splitting in lattice units, A = 5 x 10 



units this is 100(14) MeV, which is is in good agreement with the experimental 
value of 97(6) MeV. This also means that as determined from the K mass is 
appropriate to the Bs system. 



4.3.2 Hyperfine splitting 

For the hyperfine splitting we use the ratio fit method as described in section 



for the splitting between finite and zero momentum correlation functions. The 
fit interval is chosen to be tmin/tmax = 9/25. An effective mass plot is shown 
in figure |. The fit results are given in table There is no visible dependence 
on the light quark mass. This also holds for the experimental results, which are 
46.0(6) MeV for the B*-B splitting and 47.0(2.6) MeV for the B*-Bs splitting. 
We find our simulation results for the splitting to be approximately propor- 
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Figure 4: Effective mass of the ratio of the correlation functions of the B* and 
the B at arrin = 2.0, k = 0.1432. The sohd hne denotes the fit. 
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Figure 5: B* — B sphtting at ft = 0.1440 plotted against the inverse pseudoscalar 
meson mass. Squares denote NRQCD results and circles, UKQCD results using 
the clover formulation for the heavy quark, plotted against the inverse dynamical 
meson mass [|1^. The dotted lines show correlated fits of the data to linear 
functions. Errors in the meson mass are not shown for clarity. 
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tional to the heavy meson mass: 



AEocl/Mp. (29) 

This is illustrated in figure |[ Using a correlated fit to the results at k = 0.1440 
it is found that the splitting extrapolates to aAE = 0.001(2) as Mp oo. This 
is consistent with zero, as it should be. The fit result for the slope is 0.043(6). 
From a fit to the simulation results at k = 0.1432 one obtains very similar values 
for the static extrapolation and the slope. Converting the value corresponding 
to the physical B mass into physical units, we obtain 36 MeV. This is relatively 
close to the experimental value, but slightly low. This might be an effect of 
quenching, which is expected to decrease the wave function at the origin and 
thus the hyperfine splitting. The error is 5 MeV if uncertainties in the lattice 
spacing are not taken into account; if they are included one obtains an error of 9 
MeV. Note that errors in effectively appear squared in the hyperfine splitting, 
as they affect both the value of the splitting and the meson mass at which the 
splitting is taken. Previous studies by UKQCD using a relativistic quark 



action for the heavy quark have given much lower results, which are incompatible 
with experiment, when plotted against the 'static' meson mass. As we show 
in figure ||, their results, using the same ensemble of gauge field configurations 
and light propagators and clover fermions without tadpole improvement for the 
heavy quark, are compatible with ours when plotted against the kinetic mass [|l9[] . 
However there is a systematic uncertainty in the determination of the kinetic mass 
for clover fermions of about 20 % around the B, as well as statistical errors that 
are larger than 30 % in the region of the B. 



4.3.3 Ab- B and Kb - B splittings 

An example for the effective mass of the correlation function of the A;, baryon is 
shown in figure ^. The bare binding energy of the A^ is determined from a fit 
of the correlation function to a single exponential. The interval is chosen to be 
tmin/tmax = 7/13 for K = 0.1440 and tmin/tmax = 7/16 for K = 0.1432. The results 
are given in table The splitting between the A^ and the meson is determined 
by calculating the bootstrapped difference between the ground state energies of 
both states. On our ensemble there are large error bars on the baryon energy. 
From HQET we expect that the splitting between the A;, and the spin average of 
the B and the B*, B, should be fairly independent of the heavy quark mass. In 
the Ab the b quark couples to a spin zero system of light quarks. In the B and 
the B*, b couples to a light quark of spin 1/2, giving a hyperfine splitting which 
is a I/ttlq effect. This expectation is fulfilled well for our results, as is shown 
in figure ^ and is borne out by experiment, AE{Ab — B) being 330(50) MeV 
and AE{Ac — D) being 313(1) MeV. Performing a correlated fit of the splittings 
at all our heavy quark masses to a constant C, one obtains C = 0.20(3), with 
a psr degree of freedom of 0.01. The mass independence of this splitting 
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K = 0.1432 


K = 0.1440 




n E ■ 


Q 


n E ■ 


Q 


1.71 


0.82(1) 


0.25 


0.76(2) 


0.41 


2.0 


0.82(1) 


0.31 


0.76(2) 


0.43 


2.5 


0.82(2) 


0.31 


0.77(2) 


0.47 


4.0 


0.82(2) 


0.34 


0.77(2) 


0.40 



Table 10: Fit results for the ground state energy in lattice units. 





aAE{Ab - B) 


aAE{Ab - B) 


1.71 


0.21(7) 


0.19(6) 


2.0 


0.22(4) 


0.20(5) 


2.5 


0.22(4) 


0.20(4) 


4.0 


0.21(4) 


0.20(5) 



Table 11: Chirally extrapolated splitting between the A;, and the B and between 
the Ab and the spin average between the B and B* , in lattice units. 
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Figure 6: Effective mass of the A;, correlation function at awPq — 2.0, k — 0.1432. 
The solid hne denotes the fit. 
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Figure 7: The chirally extrapolated splitting between the A;, and the spin-average 
of the B and B*, plotted against the inverse spin-averaged meson mass. The 
dotted line shows a correlated fit of the data to a constant. 



makes it a good quantity to calculate and we could actually use it to extract a~^. 
Using the experimental value for the splitting at the B, we obtain = 1.7(2) 
GeV. Converting the constant C into physical units, using a lattice spacing of 
= 2.05(10) GeV, one obtains 0.41(9) GeV. Our central value is slightly higher 
than the experimental result, albeit with large errors. It is possible that the Af, 
suffers finite volume effects on these lattices. 

C is also our prediction for the static Ab — B splitting, which is in good 
agreement with the static result of the UKQCD collaboration at /3 = 6.2 of 
0.42 ^ ^g(stat) 1 l(syst) GeV 0. Our correlation functions for the static A;, 
are rather noisy at this level of statistics and it is problematic to extract baryon 
energies from them. 

Other groups extract results for the Af, — B splitting. This is expected to be 
more dependent on the heavy quark mass. Experimentally, one has AE{Ah—B) = 
0.36(5) GeV and AE{A^ - D) = 0.416(1) GeV. We work at the b quark mass 
and can extract a splitting of 0.43(10) GeV. In fact there is no noticeable mass 
dependence in our Ai, — B splitting in table [Tl|, but our error bars are larger than 
the hyperfine splitting, so at our level of statistics we could not clearly detect it. 
However, without chiral extrapolation the mass dependence is somewhat more 
obvious. Other groups have used a relativistic quark action and extrapolated 
from moderate quark masses to the b (for a compilation of results see |^). This 
is dangerous, particularly if such an extrapolation has been shown to give an 
incorrect hyperfine splitting. We believe that there could be significant systematic 
errors attached to this procedure. 
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4.3.4 P - S splitting 



The observation of orbitally excited mesons, generically called B** states, has 
been reported recently by DELPHI and OPAL ||2^, They study decay modes 
to Btc and B*n and are not able to clearly separate the expected 0"*", l"*", 1"^ and 
resonances, the first two expected to be rather broad {jg = 1/2) and the latter 
two, narrow {jg = 3/2). A mass of 5.73(25) GeV is quoted for the cross-section 
weighted mean mass of B** states by DELPHI. This gives some ambiguity when 
we wish to compare lattice calculations with experiment. The ideal quantity to 
calculate is the spin-averaged splitting, B** — B. The experiments above give a 
B** — B splitting of 419(25) MeV without spin-averaging. Thus individual B** 
states, and hence the spin-average may differ from this value by the splittings 
between P states, and it is unknown how individual states contribute to the 
experimental number. A splitting of around 50 MeV between the jg = 3/2 and 
the jg = 1/2 states has been suggested p4 
narrow 2"*" and l"*" states have been clearly seen and yield a D{1 



For the case of charm the two 
D splitting 



of 450(3) MeV. This is very similar to the number above and indicates that the 
heavy quark mass dependence of the P — S splitting isn't large and will only be 
visible once the separation of P states is clear. 
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Figure 8: Effective mass of the ^Pi correlation function, arriQ = 2.0, k = 0.1432. 
The solid line denotes the fit. 

For our P state correlation functions, we are only able to extract a good signal 
from the ^Pi. The ^Pq, ^Pi and are rather noisier and give, where visible, a 
very similar effective mass to that of the the ^Pi. We have not considered the 
^Pi/^Pi cross-correlations, so it is not clear whether the mass we are extracting 
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K = 0.1432 


K = 0.1440 




-^sim 


Q 


-^sim 


Q 


1.71 


0.70(4) 


0.39 


0.72(4) 


0.39 


2.0 


0.72(1) 


0.27 


0.72(1) 


0.25 


2.5 


0.71(1) 


0.26 


0.71(1) 


0.23 


4.0 


0.71(1) 


0.32 


0.71(2) 


0.38 



Table 12: Fit results for the ^Pi states. 





aAE 


1.71 


0.27(4) 


2.0 


0.27(2) 


2.5 


0.25(2) 


4.0 


0.25(2) 



Table 13: ^Pi - S splitting. 



from the ^Pi is that of the physical l"*" or the or some average. This the- 
oretical uncertainty, which can be improved with a higher statistics calculation 



and taking the cross-correlations into consideration |2^ , is similar to the current 
experimental uncertainty described above. 

For the heavy quark masses 2.0, 2.5 and 4.0, a fit interval oi tYnin / tmax — 3/11 
was chosen for our ^Pi correlation functions. An example for an effective mass 
plot is shown in figure ^. In the run with arriQ = 1.71 the smearing functions for 
the P states were less optimized, so there a fit interval of tmin/tmax = 8/12 was 
used. The results for the ^Pi seem to be independent of the light quark mass, 
so we choose to determine the S — P splitting by calculating the bootstrapped 
difference between the ^Pi simulation energy at k = 0.1432 and the chirally 
extrapolated spin-averaged S state simulation energy. 

In figure | we show the S — P splitting plotted against the inverse spin- 
averaged meson mass. The dependence on heavy quark mass is clearly small, 
as expected. A correlated fit to a linear function gives an extrapolation to the 
static value of 0.22(5) in lattice units and a slope of 0.1(1), compatible with zero. 
Taking our result at amq = 2.0, which is approximately at the B, and = 
2.05(10) GeV gives a splitting of 0.55(4) GeV. This is rather higher than the 
experimental value above but not inconsistent given the systematic errors arising 
from the uncertainty (both theoretical and experimental) over which B** state 
we are considering. If we use the experimental result to derive a lattice spacing, 
we obtain the rather low value 1.6(2) GeV, in fact in agreement with our result 
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Figure 9: The ^Pi — S splitting plotted against the inverse spin-averaged meson 
mass. The dotted line corresponds to a correlated fit of the lattice results to a 
linear function. 



from Af, — B. Previous results in the static approximation give a splitting 
of around 0.4 GeV for the B{1~^) — B splitting. There is some indication in our 
results that the static value might be slightly lower than that at the h mass, but 
we have not attempted to calculate this splitting in the static approximation. 



5 Conclusions 

We present the first comprehensive lattice study of the heavy-light spectrum 
at meson masses in the region of the B in the quenched approximation, using 
NRQCD and static heavy quarks. The light fermion action has been improved 
through 0{a) at tree level, and the NRQCD action has been tadpole improved. 
Work is in progress using tadpole improved light fermions with higher statis- 
tics 1^ . 



Our work has shown that with our method it is feasible to simulate heavy- 
light hadrons with a h quark directly on the lattice. Results for meson masses, 
the hyperfine splitting, Bg — Bd splitting, P wave states and heavy-light baryons 
have been obtained that are in reasonable agreement with experiment. A com- 



parison between the lattice results and experimental values is shown in figure [10 
Note that in figure |10| we have fixed the B mass to its experimental value and 
plotted splittings in physical units. This study on quenched configurations has 
been performed in parallel with an investigation with 2 flavours of dynamical 
quarks [llB, Eol , which enables us in principle to estimate the dependence of the 
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Figure 10: The B spectrum. Filled circles denote our results, where error bars 
do not take uncertainties in a"^ and the perturbative energy shifts into account. 
The dashed lines denote the upper and lower bounds on the experimental data. 
The mass of the B has been shifted upwards to match the physical value. 



results on the number of flavours. 
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